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■ ABSTRACT 

The aim of this work is to demonstrate the properties of the magnetospheric 
£ — ■ model around Kerr blackholes (BHs), so-called the fly-wheel (rotation driven) model. 

The fly-wheel engine of the BH-accretion disk system is applied to the statistics of 
QSOs/AGNs. In the model, the central BH is assumed to be formed at z ~ 10 2 
' and obtains nearly maximum but finite rotation energy (~ extreme Kerr BH) at the 

formation stage. The inherently obtained rotation energy of the Kerr BH is released 
through an magnetohydrodynamic process. This model naturally leads finite lifetime 
of AGN activity. 

Nitta et al. (1991) clarified individual evolution of Kerr BH fly-wheel engine which 

■ is parametrized by BH mass, initial Kerr parameter, magnetic field near the horizon 
| and a dimension-less small parameter. We impose a statistical model for the initial 

Q\ . mass function (IMF) of ensemble of BHs by the Press-Schechter formalism. By the 

help of additional assumptions, we can discuss the evolution of the luminosity function 
and the spatial number density of QSOs/AGNs. 

Comparing with observations, somewhat flat IMF and weak dependence of the 
magnetic field on BH mass are preferred. The result well explains decrease of very 
,H \ bright QSOs and decrease of population after z <~ 2. 

a - 
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1 INTRODUCTION that the central BH presumably gets an enormous angular 
, momentum at the formation stage of the monster. For ex- 
We discuss the evolution of QSO/AGN activities under the ample, Sasaki & Umemura (1996) showed the formation sce- 
fly-wheel (rotation driven) model which is one of the plausi- nario by using the Compton drag process as follows. After 
ble models for the powerful engine of the AGNs including a the neutralization of universe at z ~ 10 3 , density fluctua- 
rotating central blackhole (BH). This fly-wheel engine might tions grow to form the proto-galactic cloud, and each fluc- 
not be familiar comparing with the fuel engine (accretion tuation obtains the angular momentum through the tidal in- 
driven engine), however, this is very attractive because this teraction among them. At the era z ~ 10 2 , nuclear reactions 
model can explain the evolution and the lifetime of AGN occur in the central region of the rotating proto-galactic 
activities very naturally. cloud, and the matters are reionized by stellar UV radi- 
It is widely believed that recent discovery of the red tail ation. The rotating reionized matters must interact again 
of emission lines (Fe Ka) from the central region of AGNs with the uniformly distributed cosmic background radiation 
suggests that the central blackholes are quickly rotating, horn the last scattered surface. This interaction efficiently 
i.e., the monster BH should be the Kerr BH (see Tanaka extracts the angular momentum of the cloud material to the 
et al. 1995, Iwasawa et al. 1996, Dabrowski et al. 1997). radiation field by the Compton scattering, then the angular 
The innermost stable circular orbit around the Kerr BH can momentum of the matter decreases, and when it goes just 
intrude more close to the horizon than the Schwartzschild below a critical value of the angular momentum, the matter 
hole (non-rotating BH) of the same mass. Hence the line suddenly collapses and makes the quickly rotating monster 
emission from the innermost region of the accretion disk blackhole. 

should be considerably red-shifted by the gravitation, and In tMg scenari0i the initial angu i ar momentum of the 

makes the red tail. central hole shall be very close to the maximum value (say, 

From the theoretical point of view, the majority believes a ~ m where a is the Kerr parameter and m is the mass 
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of the hole) which the Kerr hole can hold. Hence it seems 
reasonable to suppose that the central blackhole is similar 
to the extreme Kerr hole at the formation stage. Such holes 
have an enormous rotation energy (~ 10 54 J ~ 10 39 Wx 10 9 yr 
for the hole with mass ~ lO 8 Af0). This is enough to explain 
the total energy release of AGNs. 

There are two different types of engines for energy pro- 
duction at the blackhole-accretion disk systems. The first 
is the well-known "fuel engine" , which is the accretion pow- 
ered engine. This is the major one which has been frequently 
adopted to explain the AGN activities. The fuel engine acts 
by a process to convert the gravitational energy released 
from the infalling matter to the radiation. The standard 
disk model (Shakura & Sunyaev 1973) is the representative 
of this model. 

The second is the "fly-wheel engine". This is the rota- 
tion powered engine. While the fly-wheel engine is not so 
familiar in the field of AGNs, this engine is as powerful as 
the fuel engine, and has very interesting features as discussed 
in the following sections. In contrast with the fuel engines, 
the energy source is the rotation of the Kerr BH itself. Of 
course, the rotation of the accretion disk also can be another 
energy source. However our scope is focused to the case in 
which the rotation of the BH is energy source, because, the 
rotation energy of the disk is supplied by accretion, so that 
there is an apprehension of confusion of two engines. The 
author strongly hopes to introduce this fascinating engine 
to researchers working in the field of AGNs. The compari- 
son of the fly-wheel engine and the fuel engine is discussed 
in section 2. 

Let us summarize the properties of the fly-wheel engine. 
The idea to extract the rotation energy of the Kerr holes is 
firstly proposed by Penrose (1969). When an incident parti- 
cle into the Kerr hole splits into two parts inside the ergo- 
sphere with a very large relative velocity, one particle can be 
thrown into the negative energy orbit falling to the hole and 
another particle escapes outward with larger energy than 
the initial energy of the incident particle. In this case, the 
rotation energy of the hole is reduced by the infall of the 
negative energy particle, and reduced energy is carried by 
the escaping particle. This is well-known "Penrose process" . 
Unfortunately it is pointed out that the Penrose process is 
not effective for astrophysical problems (see Bardeen et al. 
1972) because the critical value of the relative velocity in 
order to realize the negative energy orbit is close to one 
half of the light velocity. Such a large relative velocity can 
be realized only by nuclear reactions of particles and may 
not be achieved by usual dynamical processes, e.g., the tidal 
disruption of the accreting matter. 

The electromagnetic mechanism extracting the rotation 
energy of the Kerr BHs is firstly proposed by Blandford & 
Znajek (1977). This is well known as the "magnetic breaking 
process" or the "BZ process". In their study, the magneto- 
sphere is supposed to be "force-free" (strictly saying, this 
is "magnetically dominated"), and they clearly showed the 
energy extraction in the form of the Poynting flux when 
the rotation speed of BH is greater than it of the magneto- 
sphere. An extension of the magnetic breaking process to the 
full MHD (magnetohydrodynamic) system was performed 
by Takahashi et al. (1990) as an elementary process of the 
BH engines. By precise analysis of the MHD accretion flow 
onto the Kerr BH, they succeeded to clarify the condition 
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Figure 1. Observed luminosity function of QSOs from Boyle et 
al. (1991). QSO LF for z < 2.9 in a q = 0.5 universe. The error 
bars are based on Poisson statistics. The dotted lines indicate the 
derived model fit to the data. 

to realize the "negative energy MHD inflow". This process 
is named as the "MHD Penrose process". Nitta et al. (1991) 
studied the magnetospheric structure filled with the trans- 
magnetosonic MHD inflow onto the Kerr hole, and applied 
the MHD Penrose process to the problem of the individual 
evolution of AGNs. The result of this work is briefly reviewed 
in section 3. 

Recently, the BZ process is speculatively staged again 
as the elementary process of the 7-ray burst (GRB, see 
Paczynski 1998). In this case, the rotation energy ~ 10 47 
[J] of nearly maximum rotating Kerr BH of mass ~ 10iW© 
is considered to be extracted by very strong magnetic field 
~ 10 11 [T] in a few seconds. The extracted Poynting energy 
is expected to produce the ultra relativistic wind with the 
Lorentz factor > 10 2 . Of course, "magnetically dominated" 
assumption of original BZ process is too simple to treat the 
wind acceleration, and it should be extended to full MHD 
fly-wheel model. The fly-wheel model is still unclear, but 
it should be one of fascinating process to unify physics of 
quasars and micro-quasars. 

In this paper, the result of Nitta et al. (1991) for individ- 
ual evolution of fly-wheel engine is applied to the statistics 
of ensemble of QSOs/AGNs and compared with the observa- 
tion. Figure |l] shows observation of the luminosity function 
(LF) of QSOs. Figure shows observation of the evolution 
of the spatial number density of QSOs. Our attention will be 
focused to explain the evolution of QSOs/AGNs in a range 
< z < 5 by a mechanical process. 

In order to discuss the physical process of the plasma 
inflows and the magnetospheric structure of the Kerr BH, 
we suppose general relativistic, stationary and axisymmetric 
ideal cold MHD flows. In this case, MHD equations reduce 
to well-known basic equations: the Bernoulli equation and 
the Grad-Shafranov equation (see Takahashi et al. 1990 and 
Nitta et al. 1991) with constants of the motion. By using 
these basic equations, we discuss the properties of the fiy- 
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Figure 2. Observed evolution of QSO population from Shaver et 
al. (1996). Space densities, comoving, normalized to z = 2 ~ 3 
are plotted as a function of red shift for flat-spectrum radio-loud 
quasars with > 1.1 X 10 27 [W Hz -1 sr _1 ] (large filled circles). 
The optically selected quasar samples of Warren et al. (empty 
squares), Schmidt et al. (empty circles), Kcnnefick et al. (dotted 
square), Hawkins & Veron (small filled circles), and Irwin et al. 
(combined with Miller et al. (empty triangles)), the flat and steep 
spectrum sources studied by Dunlop & Peacock (shaded areas), 
and the flat spectrum radio-loud quasar sample of Hook et al. 
(dotted line). 



wheel engine and apply it to the evolution of ensemble of 
QSOs/AGNs in section 4. We should note again that our 
primary purpose is to demonstrate the fascinating properties 
of the fly-wheel model and not to produce a serious model 
for evolution and statistics of QSOs/AGNs. 



2 COMPARISON OF FLY-WHEEL ENGINE 
VS. FUEL ENGINE 

Here we compare the fly-wheel model with the fuel model, 
and clarify the differences between them. Our attention is 
focused to the energy source, the power output and the form 
of the energy transfer. 

The most fundamental difference is in the energy source 
of them. The energy source of the fuel engine is gravitational 
energy of infalling matter released through some dissipation 
process like a- viscosity (Shakra & Sunyaev 1973). Hence this 
engine can act while the accretion is continued. The energy 
source of the fly-wheel engine is the rotation energy of the 
central spinning BH itself which can be extracted through an 
electromagnetic process (the magnetic breaking) . We should 
note that the rotation energy of the BH is obtained at the 
formation stage, and is finite. Thus the lifetime of the fly- 
wheel engine, on the contrary, must be finite. 

The output power of the fuel engine essentially depends 
upon the mass accretion rate of the infalling matter, and 
is widely variable. Upper boundary of the power approx- 
imately corresponds to the Eddington luminosity. On the 



contrary, the output power of the fly-wheel engine is deter- 
mined by the magnetospheric equilibrium. In a typical case, 
the power ~ 10 39 [W] for the BH mass ~ 10 8 M (see the 
next section), and is enough to explain actual QSOs/AGNs. 

In the fuel engine, generated power is in a thermal form 
(e.g., the standard disk model) and is immediately converted 
to the radiation from the central region. The mechanical pro- 
cess of the fueling is so complicated. We need several models 
of the angular momentum extraction for each decade of dis- 
tance from the BH. These mechanisms must be matched 
consistently, however, this is very difficult. 

In the fly-wheel engine, the extracted rotation energy 
from the BH is once stored in the form of the Maxwell stress 
of the magnetosphere. This stress causes the global electric 
current circuit in the magnetosphere, and magnetocentrifu- 
gal force drives the plasma outflow, for example, the highly 
collimated bipolar jets in radio-loud AGNs or the equatorial 
wind in BAL QSOs. The kinetic energy of the plasma out- 
flow is finally converted to radiation at some distant region 
through some emission processes (e.g., the synchrotron radi- 
ation produced by the 1st Fermi acceleration on the shock). 
The mechanical process of the fly-wheel engine is simple. We 
can clarify it by closed discussion in the vicinity of the BH. 
In addition, we should note that the fly-wheel model can 
be a simple and unified mechanism throughout the entire 
magnetosphere in a range from AU to Kpc or Mpc. 

These properties of the fuel engine and the fly-wheel 
engine are summarized in table 1. 



3 EVOLUTION OF THE POWER OUTPUT 
FROM THE KERR BLACKHOLE 
MAGNETOSPHERE 

In the MHD scheme, the magnetosphere near the horizon 
of a BH must be filled with super-magnetosonic accretion 
flow to keep the causality. The trans-magnetosonic condition 
crucially restricts the magnetospheric structure. 

Nitta et al. (1991) studied the magnetospheric structure 
of a Kerr BH filled with trans-magnetosonic accretion flow. 
The BH magnetosphere is characterized by coexistence of 
the outgoing flow and the accreting flow. In order to real- 
ize the outgoing flow and the accreting flow@simultaneously, 
they suppose the stagnation region (source region) which is 
sustained by the magnetocentrifugal force against the grav- 
ity. This is the source of ingoing/outgoing flows. The stag- 
nation region may correspond to the pair creation region 
near the outer gap (near the surface lj — Qf where Qf 
is the angular velocity of the magnetosphere and to is the 
Lense-Thirring rotation of the inertial frame: see Hirotani 
& Okamoto 1998) or the disk halo. The accretion flow starts 
with very low poloidal velocity, and accelerates toward the 
horizon, then the flow must pass through the Alfven point 
and the fast point before reaching the horizon. 

In their result, it is clarified that the strong gravity of 
BH cause the accretion flow and amplifies the magnetic flux, 
but the total magnetic flux and the particle number flux 
Nh threading the horizon are suppressed by the rotation of 
the hole, 
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Table 1. Comparison of Fly- Wheel Model vs. Fuel Model 



Name 


Fuel 


Fly-Wheel 


Energy Source 


Gravitational energy of accreting matter 


Rotation energy of Kerr BH 


Lifetime 


Eternal 


Finite Lifetime 


Energy Transfer 


Thermarized at central region 
Radiation 


Induces the magnetospheric stress 
=*>Plasma outflow propagating toward distant region 


Applicable Region 


Central activity 


Central(~AU) to Distant(~Mpc) Region 


Properties 


Needs different fueling mechanisms 
for each decade of distance 


Closed discussion around BH is possible 
Simple & unified scheme 



where denotes the total magnetic flux of the magneto- 
sphere and u>h is the angular velocity of the Kerr blackhole 
(the Lense-Thirring angular velocity at the horizon), 



Nr 



Bo 2 



(2) 

where Bq is the magnetic field at the source region of the 
accretion flow and fi is the averaged rest mass of the particle 
of the flows. We also find that one infalling particle can 
release the rotation energy of BH of the order of its rest 
mass energy. 

According to these results we can estimate the total 
power output Lbh from the rotating BH as 

0(t - tmaxim, B , e)) 



L BH (rn,e,B ;t) = P (m,B ,e) 



yl — t/Tevo(m, Bo, e) 



(3) 



where t is the time after the birth of BH, Po = m 2 £>o/e, 
T evo — e 2 /(mBo), tmax = (1 — f?)r evo , e = mQ.F (always less 
than unity, see Nitta et al. 1991), 6 is the Heaviside func- 
tion (0(x) = for x < 0, 6{x) = 1 for x > 0). We should 
note that this formula is somewhat simplified from the orig- 
inal form (see equation [5.9] of Nitta et al. 1991). Since this 
formula is derived by estimation of order of magnitude, the 
formula contains unspecified factors of order of unity. Here 
we assume these factors to be unity for simplicity. For the 
typical case, the values of Po and r evo are given as 



1(T 



10 9 [yr 



(m 
10 8 M G 



^10 8 M G 



Bp 
1[T] 

B 
1[T] 



(4) 



(•») 



In this model, initially quickly rotating BH (ojb Of) 
spins down by the magnetic breaking process and releases 
its rotation energy. When t = t max (wh = Of), output 
power is maximum, then the engine ceases to act suddenly. 

A typical sample is shown in figure |^ where Lbh is 
denoted as a function of the time after the BH formation. 
In this scheme, the extracted energy is in the form of the 
Poynting flux, and it will be converted to the thermal energy 
of the magnetospheric plasma through some dissipative pro- 
cess or the kinetic energy of outflows by magnetocentrifugal 
drive which will act outside the source region. Let us as- 
sume here that all the extracted energy converted to the 
radiation through some unspecified mechanism, hence the 
total output power should be interpreted as the bolometric 
luminosity. We should note that during the evolution, the 
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Figure 3. Light curves of the fly-wheel engine for different BH 
mass, m = 10 6 M Q , 10 7 M Q , 10 s M Q , 10 9 Af Q , 10 10 Af Q for B = 
1[T] and e = 0.1. BHs with m ~ 1O 8 M cease to act at z ~ 3 
and important for evolutionary properties at this range. 



BH mass m is nearly constant, because the time scale of the 
mass variation is much longer than it of the angular momen- 
tum variation which determines the time scale of evolution 
(see Nitta et al. 1991). 

Here, the evolution of power output of individual Kerr 
BH magnetosphere has been clarified. We now try to ap- 
ply this result to the statistical discussion for ensemble of 
QSOs/AGNs. 

It should be noted that the strength of the BZ process 
depends crucially on that of the magnetic field threading the 
horizon. Blandford & Znajek (1977) firstly discussed the BH 
fly-wheel engine, but total magnetic flux on the horizon is 
a free parameter in their discussion based on the magneti- 
cally dominated limit. On the contrary, Nitta et al. (1991) 
using full MHD discussion enable to obtain the magnetic 
flux on the horizon as a result of the inner magnetospheric 
equilibrium based on full MHD. Instead of it, the total mag- 
netic flux of the entire magnetosphere is treated as a free 
parameter. This point will be discussed in section 6. 
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4 APPLICATION TO STATISTICS OF 
QSOS/AGNS 

General relativistic theory of QSO core formation is still an 
open question. Hence we do not have statistical properties 
of the parameters of QSO BHs. If we assume the statistical 
distributions of the BH mass m (the initial mass function), 
the Kerr parameter a (the initial Kerr parameter function) 
of seed BH and the magnetic field strength Bo at the source 
region (this should depend on the BH mass, the accretion 
rate and the dynamo theory), we can sum up the contri- 
bution of each BH over the ensemble, and can suggest the 
statistical properties of QSOs/AGNs by the Kerr BH fly- 
wheel model. 

Here we will demonstrate a preliminary application of 
Kerr BH fly-wheel model to QSO statistics. The discussion 
is based on the Press-Schechter formalism as a probable 
seed BH formation scenario. Sasaki & Umemura (1996) dis- 
cussed an additional process, the Compton drag scenario, for 
further angular momentum extraction to form the proto- 
galactic cloud to form the seed BH, and suggest the ini- 
tial mass function in figure 1 of their paper. Unfortunately, 
the distribution of the Kerr parameter and the magnetic 
field strength at the source region are not mentioned there. 
Hence we should assume the magnetic field strength Bo at 
the source region and the initial Kerr parameter a/m as 
follows. 

The magnetic filed at the source region is usually es- 
timated as Bo ~ 1[T] for m — 1O 8 M0 in order to explain 
typical QSO luminosity. This value probably depends on the 
BH mass, then we assume 



B = 1[T] x (m/10 8 M Q ) 



(6) 



We also assume the initial Kerr parameter a/m ~ 1 (nearly 
the extreme Kerr BH at the initial stage) and the small 
parameter e = mQ.F ~ 0.1. Then the power output Lbh is 
the function of m and t (see equation ^|). Thus the only we 
need is the initial mass function of BHs. 




Figure 4. Initial mass functions of BHs for n = —1.0 (dashed 
line), —0.8 (solid line), —0.5 (dot-dash line). As a result, the case 
n = —0.8 gives preferable evolution of the luminosity function 
and the spatial number density of QSOs. 



In their paper, the BH formation epoch z v ir is obtained 
from somewhat complicated procedure, however, in a BH 
mass range 1O 6 M0 < m < 10 10 Mq, z V i r varies in a very 
narrow range around z V i r ~ 200, then we neglect the mass 
dependence for simplicity and put z V i r = 200 throughout 
this paper. The resultant mass functions of BHs are shown 
in figure ^. 

We can obtain the luminosity function $ at the cosmic 
time i, 



4.1 Evolution of the luminosity function 

From Sasaki & Umemura (1996) we obtain the initial mass 
function fsn of BHs based on the standard CDM model as 



fBH{m) = 
where 



n + 3 po 
6 Mirn) 1 



-v 2 j1 



\ M c0 



(n+3)/6 



(7) 



(8) 



with M is the total (dark matter+baryon) mass of the proto- 
galactic cloud, po is the present total mass density, M c o = 
po47r(16Mpc) 3 /3. M is related with the BH mass m as 



m = rBH^ibM 



(9) 



where Qt is the fraction of the baryonic mass to the total 
mass and tbh is the ratio of the BH mass to the bakoryonic 
mass. We assume p = 6.9 x 10 10 [M o /Mpc 3 ], fl b = 0.05 
and tbh = 0.1 in this paper (we adopt a cosmological model 
with total density parameter Qo — 1 and the present Hubble 
constant Hq = 50[km/s/Mpc]). 



$(m;t) 



dn BH (m) 
dm 



dL BH (m;t) 

dL BH (m;t) 



dm 



fBH(m) 



d-L B jf (m;t) / dm 



(10) 



where Ubh (wi) is the total number of BHs having the mass 
smaller than m. Usually, evolution of the luminosity func- 
tion is parametrized by z instead of t. We adopt here the 
Einstein-de Sitter universe as the cosmological model to re- 
late t with z, 

t = t /{z + lf 2 (11) 

where to ~ 10 10 [yr] is the present time. 

Here we discuss a number of examples of the depen- 
dence of the magnetic filed Bo at the source region on BH 
mass m. The locus of the plasma source is supposed to 
several times the horizon radius (~ me -2 ' 3 ~ 4.6m for 
e = 0.1) where the pair creation seems to be effective due to 
the outer gap model (Hirotani & Okamoto 1998). The case 



Bi 



o oc m 



-1/2 



(£ = —1/2) is so-called the Eddington value. 
This formula is derived as follows. Based on the spherical ac- 
cretion with the Eddington accretion rate, we suppose the 
equipartition condition of energy density between the gravi- 



© 0000 RAS, MNRAS 000, 000-000 



6 Shin-ya Nitta 



-44 - 




35 36 37 38 39 40 41 

log(L S tf/[W]) 

Figure 5. Evolution of luminosity function of the fly-wheel engine 
for C = and n = —0.8. The brighter-end of the curve once lifts 
up for 5 > z > 2 then it drops and shrinks for 2 > z > 0. This 
tendency is consistent with observations. 



-40 



-60 - 



-80 



O 1 

1? -1 10 2 



-1.2 10 2 



-1.4 10 2 





-1 — 1 1 1 — 1 — 1 — 1 — 1—1 — 1 — 1 — 1—1— 


- 










\\ 50 




100 \ 



34 



36 



38 40 

log(W[W]) 



42 



44 



Figure 6. Evolution of luminosity function of the fly-wheel engine 
for £ = and n = —0.8 (zoom-out) in the period 100 > z > 5. 
The brighter-end of the curve considerably lifts up and shrinks 
by evolution. 



tational one and the magnetic one. However, in this case, the 
lifetime of the fly-wheel engine is determined independent of 
m from equation (^). We can easily imagine the evolution 
of the luminosity function of this case. The curve moves 
rightward holding its shape corresponding to the increase of 
output power, and arrives at the explosive stage, and then, 
all the BH engines simultaneously cease their activity. This 
is trivially inconsistent with the observation (see figure [l]), 
hence we reject this case. 

The magnetic field of the inner region of the Shakura- 
Sunyaev's accretion disk is evaluated as follows. They as- 
sumed the equipartition between the magnetic energy and 
the thermal energy. For the radiation pressure supported 
case, £ = —1/2 (similar to the Eddington value). This case 
is also unsuitable as discussed above. For the gas pressure 
supported case, £ = —1/20 assuming that the mass accre- 
tion rate is proportional to m (like the Eddington limit). In 
our discussion, since the dynamic range of the BH mass m 
is 1O 6 M0 < m < 1O 1O M0, merely in the dynamic range of 
the fourth order of magnitude, the dependence £ = —1/20 
means almost Bo ~ const. Hence we can assume that the 
mass dependence of the magnetic field Bo at the source re- 
gion is presumably very weak, —1/2 < f. In the following 
discussion, we assume £ = and set So = 1[T] independent 
of m. 

Evolution of the luminosity function is shown in figure 
|E| ^| for the typical case £ = 0, n = —0.8. In case f = 0, 
the initial luminosity is proportional to m 2 , hence the lumi- 
nosity function at the formation stage z ~ z V i r remarkably 
reflects the initial mass function (see figure ^|) . The curve of 
the luminosity function consists of monotonically decreasing 
power-law slope and exponential cut-off. 

For £ = 0, the lifetime of the fly-wheel engine is a de- 
creasing function of m (see equation []). When the BHs of 
mass, say, m = mo approach to the explosive stage t = 



tmax, the luminosity L B H(rn , t) and dL B H{m;t)/dm\m=m 
quickly increase. Until this time, the BHs having mass 
m > mo have already died. Hence the brighter-end of the 
luminosity function quickly extended in the lower-right di- 
rection (see eq. (l(|). If the initial luminosity function has 
a sufficiently steep sloop, the evolution results a lift up at 
the brighter-end. On the other hand, if the initial luminos- 
ity function has a fiat sloop, at the first, the brighter-end 
lifts up, however, the brighter-end reach the junction point 
between the power law part and the exponential cut off part 
of the luminosity function, the initial curve is bend in lower 
direction by the evolution (see figure |E| for n = —0.8). This 
behavior seems to be plausible to explain the observed evo- 
lution of the luminosity function shown in figure ^ qualita- 
tively. The criterion of the steepness is n ~ 0.8 for £ = 0. 

This evolution is just the effect of individual evolution 
of each fly-wheel engine. The time scale of the evolution is 
a decreasing function of the BH mass m, and the individ- 
ual luminosity is an increasing function of m. Hence, the 
brighter-end corresponds to massive and short lifetime BHs. 
Massive BHs quickly evolve to the explosive stage (lift up 
the curve at the exponential cut off part or more steepen at 
the power law part), and then they cease to release energy 
(shorten the curve). 

From the observational studies of the QSOs number 
counting, the bending of luminosity function curve has been 
pointed out (see, e.g., Boyle 1993, Pei 1995). The functional 
form of the curve is guessed as double power law or power 
law with exponential cut-off. Boyle (1993) argued the evolu- 
tionary motion of bending point (see figure 2 of Boyle 1993). 
The result of the fly-wheel model suggests that the bending 
point is determined by the initial mass function and does 
not move during evolution. 
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4.2 Evolution of the QSO population 

Similarly we can discuss the evolution of QSO spatial num- 
ber density (or usually called as "population") as a func- 
tion of z. At the time t, the BHs of mass m with t < 
tvir + tmaxirn) are still alive. This condition gives the up- 
per boundary m u (t) of the mass of the active BHs, because 
tmax{fn) is an increasing function of m in case Bo oc m°. 
From the condition 



tvir ~\~ (1 

we obtain 
m u (t) = 



)T RVO {m) > t 



lO 8 [M ]fe) 2 (^)" 2 (l-^)(^) 



(12) 

(13) 
(14) 



Here after we treat only the active BHs (m < m u (t)). 

Next we assume the detection limit Lum(t) of obser- 
vation on the bolometric luminosity, and let us count the 
number of BHs satisfying 



(15) 



as QSOs. Note that instead of that we cannot discuss the 
spectrum of released energy, we can treat only the total out- 
put power of the fly-wheel engine. The released energy by the 
fly-wheel engine produces, at the first, the plasma outflows, 
and finally, it is supposed to be converted to the radiation 
through some physical processes like the Fermi acceleration 
on the shock surface. Hence our attention should be focused 
to the bolometric luminosity under an assumption that re- 
leased energy is completely converted to the radiation. Of 
course, actual observational segregation of QSOs from other 
objects is based on multicolor spectroscopy, but, unfortu- 
nately we cannot argue any more than the bolometric lumi- 
nosity here. 

Models of concrete functional form of Lu m (t) will be 
given later. The condition ( |l5|) gives the lower boundary 
mi(t) of the BH mass which can be detected as QSOs at 
the time t. The relation LBH(jn;t) — Lu m (t) reduces to an 
equation for m, 

B (r 



* m 4 /e 2 + Llim {tf t B Q {m) 2 m/e 2 - L Um (tf = . (16) 



This equation can be easily solved numerically. Especially, 
for case £ = (see equation []), this equation can be reduced 
to the 4th order algebraic equation, 

fcim 4 + k-2(t)m + kz(t) 







where 
fcl = 



10>rg 2 /s 2 A4] (£f 2 (§) 



k 2 (t) 



^ 1 1 Tn 

{tftBt/e 2 



(17) 

(18) 
(19) 

(20) 



= L lim (tf ( T | F )lO- 17 [l/yr/M ](^ T )- 2 (^) 2 , (21) 

and ks(t) — —Li im (t) 2 . We can obtain unique real-positive 
root of this equation as mi(t). Let us sum up the number of 
BHs in the range mi(t) <m< m u (t), 



rm u (t) 

n Q so{t) = I f B H(m)dm. 



(22) 
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Figure 7. Evolution of QSO population for n = —0.8. QSOs 
with bolometric luminosity greater than Li im (= 10 38 ' 9 [W] for 
solid line, 10 37 ' 7 [W] for dashed line) are counted. Decrease of the 
population after the peak z ~ 2 for solid line is clearly shown as 
a result of the evolutionary scenario of Kerr BH fly-wheel engine. 
The dot-dash line denotes the population of QSOs with the energy 
flux F > 6.6 X 10 -16 [W/m 2 ]. Decrease of the population at high 
2 is more significant than above cases. 



This is the spatial number density of QSOs 
The results of simple cases Lu m (t) as 

10 38 9 [W] (solid line) 



(23) 



Llim (£) 

and 

L Hm (t) = 10 37 - 7 [W] (dashed line) (24) 

are shown in figure (?]. These values correspond to the lumi- 
nosity at absolute magnitude M = —26 (typical value for 
QSOs) and —23 (typical value for AGNs), respectively. 

Figure (?| for n = —0.8 shows very plausible evolution 
consistent with observations, in a qualitative sense, but the 
plot for n = — 1 (unplotted as figure) shows rather sudden 
decrease after the peak (z < 3). This is due to the difference 
between these cases of evolution of the luminosity function 
in this range of z. 

In the actual observation, the detection limit Lu m 
should be an increasing function of the look back time to — t 
or the red shift z, because the detection limit corresponds 
to the limit on the energy flux F like F > Fu m . We have 
also tried more realistic function of the detection limit, 



L Hm {z) = F llm 47v[VTT^ - l] 2 (l + z)/{H /cY 



(25) 



where Fn m = say 6.6 x 10 [W/m ] (apparent magnitude 
— 18.9) is the detection limit on the energy flux, 
Ho = 50[km/s/Mpc] is the Hubble constant at the present 
epoch and c is the speed of light. This is the formula of the 
translation of the energy flux Fi im to the luminosity Li im 
based on the Einstein-de Sitter universe. The result is plot- 
ted as the dot-dash line in figure Q This value of Fu m is a 
tentative and artifical one chosen to intersect the solid line 
and the dot-dash line in z = 2.5. 
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The conversion from the detection limit on the flux Lu m 
to the limit on the bolometric magnitude mi im is based on 
the relation that 

m Um = -2.5 log (26) 

where F = 2.48 x 10~ 8 [W/m 2 ] (see Allen 1973). The value 
of the bolometric magnitude mu m = 18.9 adopted here 
might be somewhat brighter than the limit in the magnitude 
of actual QSO survey with a specified band of wave length. 
However, we should note that the conversion (the bolomet- 
ric correction) between the magnitude used in observations 
(e.g., the V-band magnitude m v or the B-band magnitude 
mb) and the bolometric magnitude for AGNs is not clear. 
Thus, it is not worth to argue a detailed estimation of the 
detection limit on the bolometric magnitude here. 

At the range with large z, say z > 3, the population of 
the result plotted as the dot-dash line decreases considerably 
as z increases, comparing with the previous result (the solid 
line) in figure [?]. This is just as expected, because L\ im is 
an increasing function of z in this case. About the dot-dash 
line, in the range with small z, say z < 2, the population is 
over estimated, because Lu m which corresponds to a fixed 
apparent magnitude mu m (= 18.9 in the case in figure ^) 
is too small. Such faint nuclei should not to be classified as 
AGNs. To be more realistic, the detection limit L; im should 
be switched from (^) to (^) at the crossing point [z = 2.5 
in figure H) of these curves as z increases. 



5 SUMMARY 

The purpose of this work is not to present a realistic model 
which can precisely explain the observation, but to demon- 
strate properties of a fascinating mechanical model, i.e., the 
Kerr BH fly-wheel model, and to submit a physical scenario 
of the evolution of QSOs/AGNs not as a speculation, but as 
a result of the mechanical process. 

We have proposed a magnetohydrodynamic model for 
'engine' of QSOs/AGNs: the Kerr BH fly-wheel model (see 
section 3 of this paper and Nitta et al. 1991). This engine 
is driven by the rotation of BH. The rotation energy of BH 
is extracted by an electromagnetic process (magnetic break- 
ing). The extracted energy is once stored in the magneto- 
sphere in the form of the Maxwell stress, and then it will 
produce the plasma outflows. 

One might mislead that the fly-wheel engine is the 
mechanism only for the radio-loud activity because the re- 
leased energy produces the outflow. This might cause from 
a strong impression that highly collimated bipolar jets make 
the double radio lobes. However, from the observational 
point of view, the presence of outflows are required for other 
kind of AGNs. For example, it has been established that, 
BAL QSOs also have outflows (the "disk wind" nearly in 
the plane of the disk, see Murray et al. 1995). From the 
theoretical point of view, mechanics of the global structure 
of outflows is argued in many literature. The produced out- 
flows will show a wide variation of global structures, e.g., 
the bipolar jets of radio-loud QSOs/AGNs, the equatorial 
wind of BAL QSOs, or more (see Nitta 1994). Thus we can 
say that the radio-loud activity is simply one possibility of 
the fly-wheel engine. 



In the fly-wheel model, we can clarify the proper- 
ties and the evolution of individual engine (see section 3) 
parametrized by BH mass m, initial Kerr parameter a, mag- 
netic field Bo at the source region and a small dimension- 
less parameter e. These engines are assumed to correspond 
to QSOs/AGNs. If we obtain the statistical properties for 
these parameters, we can discuss the statistics of ensemble 
of QSOs/AGNs. Here we adopt the Press-Schechter formal- 
ism for the mass distribution. In the scenario of this work, 
Kerr BHs are supposed to form at z = z V i r = 200 with 
nearly maximum angular momentum a ~ m at the forma- 
tion epoch. The BH mass is assumed to be 10% of the total 
baryonic mass of the proto-galactic cloud. Since the mag- 
netic field Bo seems to be related with the BH mass m, we set 
Bo oc to . As a result, very weak dependence > £ > — 1/2 
is preferred for the consistency with observations. We as- 
sume C = and Bo = 1 [T] to obtain the figures. The small 
parameter e is determined by the physics of plasma injec- 
tion process, e.g., the pair plasma production or the overflow 
from the disk halo. Since we do not have widely accepted 
standard theory of it, we assume as e = 0.1 (the distance of 
the source region is several times the horizon radius). Thus 
the evolution depends only on the BH mass m. 

We have discussed the evolution of the luminosity func- 
tion and the spatial number density in a period < z < 5, 
and made a qualitative comparison with observations. In the 
typical case n = —0.8 and £ = 0, we obtain the evolution of 
the luminosity function and the spatial number density with 
a plausible behavior in the period < z < 5 consistent with 
observations. The brighter-end of the luminosity function is 
lifted up for z > 3, then it drops and the curve changes to 
be short and more steep for < z < 3 as shown in figure 
H ^. In accordance with this behavior, the spatial number 
density evolves as shown in figure [?]. 

We should note that these characteristic evolutions ob- 
tained in this paper are derived from the evolution of the 
individual magnetospheric structure in the vicinity of the 
Kerr BH. This individual evolution is not a speculation but 
the result based on the MHD picture. In the previous works, 
e.g., Pei (1995), the evolution of individual AGNs is simply 
an assumption without any mechanical scenario. We have 
tried to join the intrinsic Kerr BH magnetospheric physics, 
i.e., the fly-wheel model with observational facts, and we 
have succeed to present a mechanical model of the evolu- 
tion, at least qualitatively. 

In order to explain the observational facts of < z < 5, 
somewhat flat mass function of BHs (n ~ —0.8 in equa- 
tion and a weak dependence of the magnetic field at the 
source region on the BH mass (£ > —1/2 in equation ^|) 
are preferred. These values of n and £ should be determined 
through an extra physics, however we do not have widely 
accepted model of them, so we have treated them as free 
parameters of our picture. 



6 DISCUSSION 

6.1 Simplifications in our model 

For simplicity, we assumed the BH formation epoch as 
z V i r — 200 independent of BH mass m. From the Comp- 
ton drag model (see Sasaki & Umemura 1996), seed BH can 
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formed only at the epoch in which the background photon 
density is sufficiently high. The formation epoch should be 
2 > 10 2 . Of course, in the actual case, z V ir will depend on 
m. However, when we translate the red shift z to the cosmic 
time t by virtue of the Einstein-de Sitter model, the varia- 
tion around z ~ 200 corresponds to the order of ~ 10 6 [yr], 
and is negligible comparing with the epoch around z ~ 3 
(~ 10 9 [yr]) in which we are interested. Hence, it seems rea- 
sonable to suppose that z V i r ~ const, independent of m. 

We assumed that the dependence of the magnetic field 
i?o at the source region on the BH mass m as Bo oc m . In 
the realistic case, Bo should depend not only on m but also 
on the mass accretion rate. This problem is very difficult and 
still opened at the present time as discussed in subsection 
6.3, however, to discuss this problem as a whole is beyond 
the scope of this paper. 

We assumed that the initial Kerr parameter a as a/m ~ 
1. Bicak & Dvorak (1980) showed that extreme Kerr hole 
does not posses the magnetic field threading the horizon. 
This means that the magnetic breaking process can not ex- 
tract the rotation energy from extreme Kerr holes. Hence we 
can not set the initial Kerr parameter as a/m = 1. However 
we should note that exact value of the initial Kerr parame- 
ter is not essential. Even if a/m =, say, 0.9, 0.5 or 0.3, the 
explosive epoch t max shifted by only a factor of the order of 
unity. Such magnitude of ambiguity does not matter in our 
discussion based on an order estimation. 

The initial mass function or the initial Kerr parameter 
function of the "proto-galactic cloud" have been discussed 
in some literature (e.g., Sasaki & Umemura 1996 and Susa 
et al. 1994), but the general relativistic dynamical process 
of the BH formation from the proto-galactic cloud has not 
been solved. The initial Kerr parameter function may not 
be important comparing with the initial mass function as 
discussed in above paragraph. Hence we have supposed the 
initial mass function of the BH from it of the proto-galactic 
cloud by an assumption that 10% of baryonic mass collapses 
to form the seed BH. 

We assumed that the dimension-less small parameter 
e = mflp = 0.1 in the calculation. This value corresponds to 
the situation that the source region (the plasma injection re- 
gion or the pair creation region) is located at a radius several 
times the horizon radius (~ 4.6m for e — 0.1). This might 
be plausible for the outer gap model of the pair creation. 
The factor Qf of the definition of e roughly coincides with 
the Keplerian angular velocity at the source region. If the 
locus of the plasma source is fixed at a radius several times 
the horizon radius from the outer gap model, Qf depends 
only on m. In the evolutionary model of Nitta et al. (1991), 
m ~ const, during the characteristic time scale of the angu- 
lar momentum extraction which is defined as the time scale 
of the evolution. Hence we can assume that Of ~ const., 
thus e ~ const, independent of the time. 

In our model, each BH suddenly ceases to release energy 
at t — t max corresponding to the situation uih = If Of 
of a BH magnetosphere distributes in a range, say, Q,i < 
Of < ^2 as a function of the magnetic flux function, the 
fly-wheel activity gradually ceases on a field line in order 
of decreasing value of Of- If this variety of Qf does not 
result the vary e in order of magnitude, it does not affect 
our qualitative discussion. 

In order to relate the time after the formation of black- 



holes and the redshift, we need a cosmological model. In 
the main discussion of this paper, we adopt the Einstein- 
de Sitter model (see eq. [H]) for simplicity. Of course, the 
Einstein-de Sitter model is very classical, and recent obser- 
vational studies of cosmology support the model with the 
cosmological constant, i.e., the Lemaitre model. Here we 
have to estimate the difference of results between the cases 
with the Einstein-de Sitter model and the Lemaitre model. 
We provide the look-back times ti(z) and ti(z) of the epoch 
with the redshift z in the Einstein-de Sitter model and the 
Lemaitre model, respectively. If we choose the parameters 
flo = 0.1 and Ao = 0.9 in the Lemaitre model, the look-back 
time of the blackhole formation epoch (z = 10 2 ) can be es- 
timated as ti — 1.92ii. At the epoch z — 5 in which we 
are interested here, t2 = 1.83ti. These differences are sim- 
ply in the factor of the order of unity. Hence the difference 
of the results between these two cosmological models does 
not matter in our discussion because the evolutionary pro- 
cess of our discussion (see section 3) is based on an order 
estimation. 

6.2 2-types of BH engines 

In this paper, our attention is focused to the fly-wheel model. 
However, BH-accretion disk systems also seem to include 
another type of engine: the fuel (accretion powered) en- 
gine. The elementary process of the fuel engine is release 
of the gravitational energy of the accreting matter, so that 
the activity strongly depends on the mass accretion rate. 
It is widely believed that the accretion rate is roughly de- 
termined by the Eddington limit. This idea is based on a 
speculation that the regulating stage of the entire accretion 
process is the final stage, i.e., the accretion onto the BH. 
However, we do not have any theoretical conviction that 
how the entire system (i.e., a galaxy) determines the accre- 
tion rate. In another word, how the system can remove the 
angular momentum of the accreting matter to realize such 
accretion rate. 

In order to determine the activity of the fuel engine, we 
must solve the extraction of the angular momentum in a very 
wide spatial range. If we want to obtain the accretion rate of 
a range from an angular momentum extraction mechanism, 
we need the accretion rate outside this region where another 
mechanism may regulate the accretion rate, and so on. This 
endless chain seems to be hopeless to solve completely. Thus, 
it is so difficult to assemble the angular momentum extrac- 
tion mechanisms of each range into a consistent theory, that 
we do not have another way except to treat the final accre- 
tion rate onto the BH as a free parameter. This is the point 
of difficulty to make evolutionary scenario of QSOs/AGNs 
based on the fuel model. 

In the model adopted in this paper, we do not consider 
the activity of the fuel engine at all. However, the author 
believes that coexistence of these 2 types of engines is un- 
doubted. Even after the fly-wheel engine ceases to release 
energy, the fuel engine continues to act while the accretion 
is continued. If we make an appropriate assumption to esti- 
mate the luminosity due to the fuel engine, we should add 
the contribution of the fuel engine. For simplicity, if we as- 
sume a constant energy release of the fuel engine, i.e., the 
Eddington luminosity (depends only on m), this luminos- 
ity is comparable with the initial luminosity of the fly-wheel 
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engine. In this model, the fuel engine works constantly and 
the fly-wheel engine works in a period t max after the BH for- 
mation. Even after the fly-wheel engine dies, the luminosity 
does not vanish, but decreases to the Eddington luminosity. 

In section 2, the fly-wheel model is characterized as that 
the mechanism of it can be clarified by the closed discussion 
in the vicinity of the BH. However this statement might be 
somewhat exaggerate. The fly-wheel engine seems to relate 
with the fuel model at the point that the magnetic field Bo 
at the source region may depend on the accretion rate. Bo is 
provided as the magnetic field strength averaged in a larger 
macroscopic scale than the scale of turbulence generated in 
the accretion disk. Such large scale magnetic field should 
be amplified by the accretion. The accretion plasma carries 
the frozen magnetic field into the inner magnetosphere and 
compresses it. Against this process, small but finite resistiv- 
ity dissipates the magnetic field. Then the saturated level of 
the magnetic field strength is determined by the equilibrium 
of the compression and the dissipation. Unfortunately, this 
problem has not been solved yet as in the next subsection. 
Hence we have assumed that Bo depends on the BH mass m 
like equation (^|) because the accretion rate seems to depend 
on m. 

6.3 Ambiguity of the fly-wheel power estimation 

Ambiguities are still remained on the estimation of the 
power of the fly-wheel activity mainly due to the following 
two reasons. The first is due to some theoretical ambigu- 
ities of the estimation of the magnetic field strength near 
the BH. The second is ambiguity of the innermost magne- 
tospheric structure of the BH-accretion disk systems, espe- 
cially whether the field lines threading the horizon (or the 
innermost region of the accretion disk) are open toward the 
infinity or closed to be loop-like one. 

The poloidal magnetic field strength is traditionally es- 
timated by an intuition of a principle of equipartition be- 
tween the magnetic energy density and the gravitational one 
or the thermal one (see, e.g., Shakra & Sunyaev 1973). Re- 
cent numerical studies based on nonlinear evolutionary pro- 
cess of the resistive MHD seem to support the result from 
the equipartition. For example, Matsumoto et al. (1997) con- 
clude that predominantly toroidal magnetic field is amplified 
by a differential rotation of the disk and the plasma /3-value 
of p ~ 10 can be achieved. If there is significant magnitude 
of the poloidal magnetic field, the saturation level will be 
more large and j3 ~ 1 might be achieved. 

While this conjecture of equipartition is now widely ac- 
cepted, there is still room for a disagreement about this 
point. Recently Livio et al. (1999) critically assess the ef- 
ficiency of the Blandford-Znajek (BZ) process (the magneti- 
cally dominated case of the BH fly-wheel model) comparing 
with other disk activities. They reconsider the field strength 
in the innermost region of BH-disk system. In their result, 
power of the BZ process is dominated by the fly-wheel or 
the fuel (viscous heating) power of the innermost region of 
the accretion disk. The problem of the saturated strength of 
the poloidal magnetic field is still an open question. 

However we should note that the energetics of the fly- 
wheel process depend not only on the strength of the mag- 
netic field but also on the inner magnetospheric structure 
of the BH-disk systems. The extracted Poynting energy flux 



due to the fly-wheel process is carried along the poloidal 
magnetic field lines, and will be converted to the kinetic 
energy of plasma outflow at some distant region from the 
horizon. Hence only open magnetic field lines can take place 
to the energy extraction toward very distant region. For ex- 
ample, Nitta et al. (1991) give schematic figure for the inner- 
most magnetospheric structure (see figure 3 of that paper). 
In their result, the magnetic field lines connecting to the in- 
nermost region of the disk are closed (a loop-like structure 
connecting the BH and the disk), and do not contribute to 
the energy extraction. Open field lines are emanated from 
a high latitude region of the BH and the outer part of the 
disk. In this case, the discussion of Livio et al. (1999) should 
be altered. 

Thus efficiency of the fly-wheel process is closely com- 
bined with the disk dynamo process and the magnetospheric 
structure of the innermost region. These are very important 
but still open questions in the current state, and to argue 
this point would carry us too far away from the purpose of 
this paper. 

6.4 Discrimination of QSOs 

As discussed in section 4, in the fly-wheel model, we can only 
estimate the total output power of the engine, and we cannot 
discuss the spectrum of resultant radiation. Hence the only 
way to distinguish QSOs/AGNs from normal galaxies is set- 
ting a criterion, say Lum, on the bolometric luminosity. Here 
we suppose that entire released energy is perfectly converted 
to radiation, and suppose that the engines having the lumi- 
nosity greater than Lu m can be treated as QSOs/AGNs. 
This simplified procedure is obviously far from the actual 
QSO number counting studies. In future study, the problem 
of the resultant spectrum of the radiation should be solved. 
This is possible only if we solve the physics of plasma out- 
flows being generated by the fly-wheel engine. This is, need- 
less to say, one of the most difficult open questions in the 
magnetospheric astrophysics. 

From figure |^, the locus of the peak of population 
strongly depends on the criterion Lum - If we set the smaller 
Lum, the peak shifts to the smaller z. This means that if 
we survey AGNs in more deep, we will find more and more 
faint AGNs including low mass BH. These may correspond 
to Seyferts. However we should note if the mass of central 
BH is too small, the nucleus activity is dominated by its host 
galaxy. Such objects may not be classified as QSOs/AGNs. 
In this meaning, the criterion L iim = in 389 or 10 37 7 [W] 
adopted in this paper might be plausible, because these 
values dominate the typical luminosity of normal galaxies 
10 37 [W] (Andromeda galaxy) . 

6.5 Similarity of radio- loud and radio-quiet AGNs 

Similarity among all kinds of AGNs is widely accepted from 
the observational point of view. The spectral energy dis- 
tribution (SED) of radio-loud AGNs and radio-quiet AGNs 
are quit similar except the radio range (see, e.g., Elvis et 
al. 1994). We also cannot find any intrinsic difference in the 
evolution of the spatial number density of the optically se- 
lected QSOs, the flat-and-steep spectrum sources and radio- 
loud QSOs (see Shaver et al. 1996). These observational ev- 
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idences of similarity implicitly suggest the universality of 
physical process of QSOs/AGNs. 

As mentioned in section 5, the author believes that the 
fly-wheel model is applicable not only for radio-loud AGNs. 
The variation of AGNs might be caused from the variation of 
proper parameters (BH mass, BH angular velocity, magnetic 
field strength, etc.). The difference of parameters will lead 
to different structure of outflows (see Nitta 1994). We may 
expect that outflows having different structure will produce 
different types of spectrum of the radiation. However, the 
correspondence of a global structure of the outflow and a 
resultant spectrum of radiation is still unclear. This should 
be a future problem. 

6.6 Other stories of the fly-wheel activity 

There are other stories to make Kerr BH as central BH of 
QSOs/AGNs. Here we mention two of them: the merging of 
BHs and the spins up by accretion. 

Wilson & Colbert (1995) discussed the formation of 
AGN BHs by merging process. They tried to explain the dif- 
ference between radio-loud and radio-quiet AGNs. As well 
known, number fraction of radio-loud AGNs is only 10% of 
total AGNs, and radio-loud galaxies are mainly observed as 
elliptical galaxies. From these points, they supposed that 
radio-loud AGNs are merger events of BHs. Merging of 2 
galactic nuclei produce a quickly rotating Kerr BH, and 
a period after merging, the Blandford-Znajek process acts 
and shows radio-loud activity. Similarly, Moderski & Sikora 
(1996) supposed to make quickly rotating BH by very large 
mass accretion. 

The scenarios discussed in these papers are alternative 
one. There is room for another possibility that the fraction 
of radio-loud to radio-quiet may be related to probability to 
make jet-like outflow. The fly-wheel engine can form various 
structures of outflows (see Nitta 1994). If well-collimated jet- 
like structure which is nearly perpendicular to the galactic 
disk is formed, this will be observed as radio-loud one (or 
the radio galaxy, see Urry & Padovani 1995), because the 
terminal shock in the jet will locates far from the galactic 
disk due to very low ambient matter density in this direction. 
If equatorial wind is formed, this will be BAL QSO (see, e.g., 
Cohen et al. 1995 or Murray et al. 1995 ). Of course, these are 
simply speculations, because we do not have widely accepted 
physics for the structure formation of plasma outflows yet. 
Anyway, since we do not have authorized theory, we must 
try to test various possibilities. 

In some literature, the fuel process (accretion from disk) 
and the fly-wheel process (Blandford-Znajek process) are 
simultaneously considered (see, for example, Moderski & 
Sikora 1996 and Ghosh & Abramowicz 1997). In these mod- 
els, the accretion contributes to spin up the BH, but the BZ 
process suppresses it. On the contrary, Nitta et al. (1991) 
imply that MHD accretion onto the Kerr BH extracts the 
angular momentum and spins down the rotation of the Kerr 
BH when Of < uh (see section 3). This is a natural result 
of MHD accretion onto the Kerr BH (see Takahashi et al. 
1990). 

One might think this is contrary each other. However 
the author does not think so. Moderski & Sikora (1996) and 
Ghosh & Abramowicz (1997) suppose to start with slowly 
rotating BH (Of > uh), on the contrary, Nitta et al. (1991) 



suppose quickly rotating BH (Qf <S uJh)- This difference 
causes from the difference of concepts of the models. How- 
ever, we should note that this difference is essential. 

In our model, the energy source of the Kerr BH fly- 
wheel engine is inherently obtained rotation energy of the 
central BH. The origin of this energy is the tidal interaction 
during the collapse of the proto-galactic cloud with other 
density fluctuations. Once the fly-wheel engine starts to act, 
the rotation energy monotonically decreases, and stops at 
the state Q.f — uJh- 

On the other hand, in other models, accretion energy 
converts to the fly-wheel type activity, and they relate to the 
radio-loud activity. For example, central BH is normally in a 
state of slow rotation, however, if coalescence of BHs (Wilson 
& Colbert 1995) or very large mass accretion (Moderski & 
Sikora 1996) occurs, central BH spins up and the fly-wheel 
engine starts to wok. In these models, the energy source is, 
consequently, the accretion energy. 



6.7 Dormant quasar: Fornax A 

Let us notice a splendid example Fornax A (NGC1316, see 
Iyomoto et al. 1998) which seems clearly show properties of 
the fly-wheel engine. In general, we cannot observe evolu- 
tion of a galaxy because of very long lifetime of it, however 
Fornax A is a particular case which we can obtain the ev- 
idence of the evolution for recent 0.1 Gyr. This is a radio 
galaxy with double radio lobes. The nucleus should be active 
(> 4 x 10 34 [W] in 2-10 keV X-ray luminosity) at least 0.1 
Gyr ago, while the present activity is 'dormant' (2 x 10 33 [W] 
in 2-10 keV X-ray luminosity). 

We can guess the reason of it based on the fly-wheel 
model as follows. The fly-wheel engine was still active 0.1 
Gyr ago, and the nucleus ejected plasma outflows (bipolar 
jets) and made the radio lobes. At an epoch within a past 
0.1 Gyr, the fly-wheel engine ceased to work and nucleus 
became to be dormant. The radio lobe can emit radiation 
within a period determined by the Synchrotron cooling time 
without the energy supply by the fly-wheel engine. In this 
sense, the fly-wheel engine of Fornax A is not 'dormant' 
one but 'dead' one if without some mechanisms to spin up 
the central BH again. However the fuel engine can act after 
the fly-wheel activity cease. This corresponds to the present 
nucleus activity. In this case, the peak fly-wheel activity is 
an order of magnitude greater than the fuel activity. This is 
just a result of our model with e ~ 0.1. 
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